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Data are presented on the reaction J/ip — » ppr\ using 7.8 x 10 6 J/V> trig gers col- 
lected by the BEjing Spectrometer (BES). A partial wave analysis is performed. 
A clear enhancement near the pr\(j>r\) threshold is observed. It is fitted with a 
J p = i" resonance with mass M = 1 530 ± 10 MeV and width T = 95 ± 25 MeV. 
In addition, there is a peak around 1650 MeV with J p = i preferred also, fitted 
with M = 1647 ± 20 MeV and V = 145+1? MeV. These two N* resonances are 
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believed to be the two well established states, Sii(1535) and 5*11(1650), respec- 
tively. It is the first partial wave study of the production of these resonances 
from J/ip decays. 

PACS: 13.25. Gv; 14.20Gk; 13.65.+i 



I. INTRODUCTION 

Nucleons are the most common form of hadronic matter on the earth and probably in the 
whole universe. To understand the internal quark-gluon structure of nucleon and its excited 
states iV*'s is one of the most important tasks in nowadays particle and nuclear physics. The 
main source of information for the nucleon internal structure is their mass spectrum, various 
production and decay rates. Our present knowledge on this aspect came almost entirely from the 
old generation of irN experiments of more than twenty years ago [Q . Considering its importance 
for the understanding of the nonperturbative QCD Q , a new generation of experiments on N* 
physics with electromagnetic probes (real photon and space-like virtual photon) has recently 
been started at new facilities such as CEBAF at JLAB, ELSA at Bonn, GRAAL at Grenoble. 

The J/tp experiment at the Beijing Electron-Positron Collider(BEPC) has long been known 
as the best place for looking for glueballs. But in fact it is also an excellent place for studying 
N* resonances ||, especially in the mass range 1-2 GeV. The corresponding Feynman graphs 
for the N* and N* production are shown in Fig. [l]. These graphs are almost identical to those 
describing the N* electro-production process if the direction of the time axis is rotated by 90°. 
The only difference is that the virtual photon here is time-like instead of space-like and couples 
to NN* through a real vector charmonium meson J/ip. Therefore almost all channels of N* 
decays studied at the CEBAF and other r )p{ep) facilities can also be studied here. 




FIG. 1. Feynman graphs for N* and N* production from e + e collision through J/tp meson 

Among many interesting channels, the J/ip — > ppr\ is a relatively simple one to begin with. 
According to the information |I]j4|-^| from irN — > rjN and 7A — > rjN experiments, as well 
as some early quark shell model calculation 0, only iV*(1535)S , n has a big decay branching 
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ratio to Nrj, while N* (1650)Sn may have a branching ratio to Nr] of up to 10% and other N* 
resonances below 2.0 GeV have much smaller branching ratios to the Nr). 

In this paper, we present BES data and a partial wave analysis on the J/ip — ► pprj reaction. 
The iV*(1535) and TV* (1650) are observed. This is the first partial wave study of the production 
of these iV* resonances from the J/ip hadronic decays in the world. The new information on 
J/ipNN* couplings provides a new source for studying baryon structure H . 



II. BES DETECTOR 

The analysis in this paper uses 7.8 x 10 6 J/ip triggers collected by the Beijing Spectrom- 
eter(BES). BES is a conventional solenoidal magnet detector that is described in detail in 
Ref. |^]. A four-layer central drift chamber(CDC) surrounding the beampipe provides trigger 
information. A forty-layer cylindrical main drift chamber(MDC), located radially outside the 
CDC, provides trajectory and energy loss (dE/dX) information for charged tracks over 85% of 
the total solid angle. The momentum resolution is crp/P = 0.017\/l + P 2 (p in GeV/c), and 
the dE/dX resolution for hadron tracks is 11%. An array of 48 scintillation counters(with 
inner radius of 1.157 meter and outer radius 1.207 meter)surrounding the MDC measures the 
time-of- flight (TOF) of charged tracks with a resolution of 450 ps for hadrons. Radially outside 
of TOF system is a 12 radiation length thick, lead-gas barrel shower counter(BSC) operating 
in the limited streamer mode. This device covers 80% of the total solid angle and measures 
the energies of electrons and photons with an energy resolution of &e/E = 22% /VE (E GeV). 
Outside the BSC is a solenoid, which provides a 0.4 Tesla magnetic field over the tracking 
volume. An iron flux return is instrumented with three double layers of counters that identify 
muons of momentum greater than 0.5 GeV/C. 



III. EVENT SELECTION 

The T) is detected here in its 77 decay mode. Therefore, much effort has been devoted to the 
selection of the events in the 2^fpp final state. Each candidate event is required to have two 
oppositely signed charged tracks with a good helix fit in the polar angle range —0.8 < cos 9 < 0.8 
in MDC and at least 2 reconstructed 7's in BSC. A vertex is required within an interaction 
region ±15 cm longitudinally and 2 cm radially. A minimum energy cut of 60 MeV is imposed 
on the photons. Showers associated with charged tracks are also removed. 

After above selection, we use TOF information to identify the pp pairs, and at least one track 
with unambiguous TOF information is required. The open angle of two charged tracks smaller 
than 175° is required in order to remove back to back events; to remove radiative Bhabha 
events, we require {E+/P+ - l) 2 + (E-/P- - l) 2 > 0.3, where E + , P + (£L, P_) are the 
energy deposited in BSC and momentum of positron (electron) respectively. Events are fitted 
kinematically to the 4C hypotheses J/ip —> 2jpp. If the number of the selected photons is larger 
than two, the fit is repeated using all permutations of the photons. For events with a good fit, 
the two photon combination with the largest probability is selected. Fig. || shows the invariant 
mass spectrum of two gammas, we can see clear 7r° and 77 signals. Meanwhile, the events are 
also fitted to J '/tp — > jpp and 4jpp. We require 

Prob(x 2 (2 -ypp) > 4C) > Prob( X 2 hpP) , 4£7), Prob{ X 2 [2lpP) , 4C) > Prob( X 2 (4lpP) , AC) 
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to reject the jpp and ppTr ir° backgrounds. In order to suppress further the backgrounds with a 
7T°, a 5C fit for the ppr\ with 77 — > 77 is performed on the selected events. Prob(Xpp v , 5C) > 1% is 
required. This 5C fit helps improving the mass resolution for combinations of charged particles. 
There are 765 events which survive our selections. 




Myy mass(GeV/c 2 ) 



FIG. 2. 77 invariant mass spectrum after 4C fit for J/ip — > PP77 
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FIG. 3. Dalitz plot for J/tp ^ pprj (M% v v M? v ) 
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FIG. 4. prj(prf) invariant mass spectrum for J/ip — » pp?7, crosses are data and histogram the fit 



Fig. |^ is the Dalitz plot for the decay J/ip — > ppry. The corresponding prj(prf) invariant mass 
spectrum is shown in Fig. ||. There is a clear enhancement near the pr)(pr]) threshold. There 
are also some bumps around 1.65 GeV and 1.8 GeV. 



IV. SIMULATION OF SIGNALS AND BACKGROUNDS 

In order to estimate the selection efficiency, a phase space generator is used to produce a 
sample of ppr\ events with full detector simulation, it is 18% after the Monte Carlo data go 
through the same cuts as real data. 

Since 2"fpp is the final state of the reaction channel, the possible background may come from 
J/ip — > ppir , ppir Q ir°, ppw, tt + tt^tt° and K + K~ir° decay modes. For each channel, 10,000 MC 
events are generated. These MC events go through the same analysis program as for the real 
data. We find that only ppir°ir° channel gives a significant contribution to the background. The 
estimated background is 8%, 

A correct estimation of the background invariant mass shape is essential to the partial wave 
analysis. Fortunately the invariant mass spectrum of these events is of the same shape as the 
J/ip — > ppr\ phase space distribution. Fig. |5| shows the mass distributions of the background from 
J/ip — ► ppir°ir°, compared with the phase space distribution for the J/ip — ► ppij process from 
Monte Carlo data. A side-band method is also used to check the shape of possible backgrounds 
from real data; the same result is obtained. 
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FIG. 5. The pr\ invariant mass distribution of the background from J/ip — » pp-K°-K° channel (crosses) 
compared with the phase-space distribution for the J/ip — > ppq reaction (histograms) 



V. AMPLITUDE ANALYSIS 

Based on the study of pp and pr\ invariant mass distributions in our data, we are mainly 
interested in the structures at 1535 and 1650 MeV of the prj invariant mass. According to 
PDG 0, below 1900 MeV there are eight N* resonances observed by other experiments, i.e., 
Pn(1440), Sn(1535), 5 n (1650), £> 15 (1675), F 15 (1680), r> 13 (1700), P u (1710) and Pi 3 (1720) 
states, among which only S'ii(1535) and Sii(1650) were observed in the rjN decay mode. Al- 
though the two prj structures at 1535 and 1650 MeV in our data are most probably due to 
Sii(1535) and 5*11(1650) states, we perform the partial wave analysis by allowing all possible 
quantum numbers of J p = * , | and | . 

The background from multi-7r° is ~ 8% in the 5C fit. We have included a phase space 
background in the PWA fit to allow for this. 

In this analysis, we use the effective Lagrangian approach for the partial wave analysis. The 
relevant spin-1/2 interaction Lagrangians are [pd] , p^l 

Kpb = -ignPRPTRr] + H.c, (1) 

£ %* = M^^<rpr+H.c., ( 2) 

4ph = -9v R RT,PV + H.c. (3) 

where R is the generic notation for the resonance with mass Mr, P for proton with mass 
Mp and ip for J/ip with four-momentum q. The vertex coupling constants g v PR, gr R and gv R 
are parameters to be determined by fitting experimental data. Note here from the single pprj 
channel, we can only determine products of coupling constants, i.e., gr}PRgT R and g-qPRgvn- 
The operator structures for the T, T ^ and T ^ v are 

r = i, r M = 7 5 7 M , = 750-^, (4) 
T = 75, r M =7 M , T lxv =a llv , (5) 
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where ^ and (||) correspond to nucleon resonances of J p = | and i + , respectively. The 
interaction Lagrangians for J p = | and § TV* resonances are constructed similarly [131. 
The amplitudes in the PWA analysis are constructed from these Lorentz-invariant interactions 
and the iV* propagators for J/ip initial states with helicity ±1. The relative magnitudes 
and phases of the amplitudes are determined by a maximum likelihood fit to the data. The 
BW parameters for different states are fitted. The pij mass projection fitted to the real data is 
shown in Fig. ^. Fig. |^ shows the data and fit for the polar angle of the proton measured with 
respect to the beam direction. Fig. fj] shows the invariant mass spectra for pp. We now discuss 
the features of the data and the outcome of fits. 
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FIG. 6. The polar angle of proton measured with respect to the beam direction, crosses are data 
and histograms the fit 
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FIG. 7. pp invariant mass spectra for J/ip — > ppr], crosses are data and histograms the fit 
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511(1535) 

The peak around 1535 MeV near the pn threshold optimizes at M = 1530 ± 10 MeV and 
r = 95 ± 15 MeV. The data favour J p = | over others. A fit with J p = i instead gives 
In L worse by 7.5 than for ^ assignment (Our definition of In L is such that it increases by 0.5 
for a one standard deviation change in one parameter). Fits with other quantum numbers are 
much worse, e.g., lni worse by 26 for J p = | . With our 4 fitted parameters, the statistical 
significance of the peak is above 6. Oct. It is obviously the 511 V*(1535) resonance. It makes 
the largest contribution (56 ± 15)% to the ppr\ final states, the errors here and later include 
both statistics and systematic errors from the fit. Our results for iV*(1535) are consistent with 
the resonance parameters of PDG jlj and a recent detailed analysis of nN 5n partial wave by 
Vrana, Dytman and Lee [fl5"|| . 
5n(1650) 

The second peak around 1650 MeV is also fitted with a J p = \ resonance iV*(1650). It 
optimizes at M = 1647 ± 20 MeV, T = 145±fj? MeV, and contributes (24±f 5 )% to the pprj 
final states. We have tried fits to this peak with other quantum numbers. The log likelihood 
is worse by 3, 4, 5, 6, 12 for | , | , | , \ and § , respectively. Note there are two more 
free parameters for J p = ^ and J p = ^ modes than for J p = ~ modes. If we assume no 
contribution from resonance around 1650 MeV, then the mass and width of 5n(1535) optimize 
around 1570 MeV and 270 MeV, respectively, which are out of the range of PDG values. The 
log likelihood is worse by 14. 

A small improvement to the fit is given by including a J p = | resonance, which optimizes 
at M = 1800 ± 40 MeV and T = 165± *f MeV. It contributes (12 ± 7)% to the ppn final states. 
The statistical significance of the peak is only 2. Oct. We have tried other quantum numbers. All 
of them give equally good fits within one standard deviation. There is a theoretical prediction 
that the lowest-lying hybrid states should be around this energy region. 
We also include a small contribution (8 ± 4)% from the tail of Pn(1440) resonance with its 
mass and width fixed to its PDG central values for its pole position [|). It is also only a 2. Oct 
effect. 

Fig. H shows the different component contribution. 
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(2) 

An interesting result is that the C]p PR term given by Eq. (3) makes insignificant contribution 
for both Sii(1535) and 5*11 (1650). If we drop this kind of couplings for both resonances, 
the likelihood value for the fit is only worse by 0.9 for 4 less free parameters. This kind of 
couplings should vanish for the real photon coupling to NN* due to the requirement of the 
gauge invariance. Why it also vanishes for the ipNN* coupling needs to be understood. A 
theoretical calculation [|l6| assuming pure C^ PR coupling without C^ PR coupling failed to 
reproduce the basic feature of the J/ip — > pprj data. This is consistent with our observation 
that the C { ^ PR coupling dominates for both iV*(1535) and 7V*(1650). 

In the Vector Meson Dominance (VMD) picture, the virtual photon couples to the NN* 
through vector mesons, and the electro-magnetic NN* transition form factors g 7 * nn* can 
be expressed in terms of photon-meson coupling strengths C 7 y and meson- NN* vertex form 
factors g VNN , : 



g^NN*(q 2 ) = ' 



r9v- N N 
j 



(<7 2 ) 



(6) 



with 



C 7 y — 



sr 



amy 



(7) 



At q 2 = ML the J/ip meson dominates. The terms from other vector mesons are negligible. 
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From our PWA results here and other relevant information from PDG we can deduce the 
transition form factor for the time-like virtual photon to PiV*(1535) as 



\g 7 * pN *(q 2 =M|)| = 2.3 ±0.5, 



(8) 




(9) 



This is the first measurement of the form factor g~ t *pN* for a- time-like virtual photon. It 
provides a new challenge for various theoretical models to reproduce it, complementary to the 
information for a real or space- like virtual photon [17|. 



In summary, the J/ip decay at BEPC provides a new excellent laboratory for studying the N* 
resonances. We have studied the J/ip — ► pprj decay channel, and a PWA analysis is performed 



MeV with r = 95 ± 25 MeV near the rjN threshold. In addition, there is an obvious resonance 
around 1650 MeV with J p = |" preferred, M = 1647 ± 20 MeV and T = U5±f 5 MeV. In the 
higher pri(prj) mass region, there is an evidence for a structure around 1800 MeV; with present 
statistics we cannot determine its quantum numbers. 

All above analysis is the first step for us to probe N* baryons at BES. We will perform a 
detail study of N* baryons in the following J/ip decay channels: J/ip — > ppTr a , pnir~ , ppir°ir°, 
pp7r + 7r~, pprj' ', ppuj and so on. 

With the forthcoming 50 million more J/ip events in near future, more precise partial wave 
analyses can be carried out on many channels involving N* resonances and should offer some 
determinations of N* properties. A systematic experimental study of the N* production from 
the J I ip decays is underway and will provide a new domain to explore the internal quark-gluon 
structure of these excited nucleons. 
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VI. CONCLUSION 



on the data. There is a definite requirement for a J p = ^ 



component at M = 1530 ± 10 
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